Anomalous bond stretching phonons as a probe of charge fluctuations 

in perovskites 

S. Cojocaru a ' c >* R. Citro a M. Marinaro a ' b 

a Dipartimento di Fisica "E. R. Caianiello" , Universitd degli Studi di Salerno and CNISM, Unitd di ricerca di Salerno, Via 

S. Allende, 84081 Baronissi (SA), Italy 
b I.I.A.S.S., Via G. Pellegrino, n. 19 84019 Vietri sul Mare (SA), Italy 
c Institute of Applied Physics, Chisinau 2028, Moldova 



Abstract 

Important information on momentum resolved low energy charge response can be extracted from anomalous properties 
of bond stretching in plane phonons observed in inelastic neutron and X-ray scattering in cuprates and some other 
perovskites. We discuss a semiphenomenological model based on coupling of phonons to a single charge mode. The 
phonon dispersion and linewidth allow to locate the energy of the charge excitation in the mid infrared part of 
the spectrum and to determine some of its characteristics. New experiments on oxygen isotope substitution could 
allow to achieve a more detailed description. Corresponding relations following from the model can be used for the 
interpretation of experiments and as test of the model. 



1. Introduction 

A direct coupling of the bond stretching (BS) in- 
plane (C11O2) phonon modes to density fluctuations 
of strongly correlated electrons has been discussed 
as a possible reason of their anomalous properties 
[1],[2]. Doping dependent anomalies of these high 
energy (70 — 90 meV) phonon branches are observed 
in a number of perovskite materials like cuprates, 
bismuthates, manganites, nickelates [3]- [7]. In par- 
ticular, presence of a strong (more than 20%) soft- 
ening of the BS phonon dispersion towards the Bril- 
louin zone boundary in the (100) and (010) direc- 
tions contrasts a very week dispersion in the diago- 
nal direction (110) . This is paralleled by the phonon 
linewidth that reaches unusually large values too, 
10 meV, in the (100) and (010). In some materi- 
als softening can also be non-monotonous and one 
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observes a minimum at intermediate momenta, e.g. 
q ~ 0.3 r.l.u. A number of theories and models 
have been developed to explain the properties of 
BS phonons. To mention a few possible scenarios, 
such as incommensurate charge density wave insta- 
bility or stripes, acoustic and optical plasmons or 
unusual dielectric properties (overscreening, when 
the dielectric constant becomes negative at some 
momenta) as well as calculations based on ab-initio 
LDA, e.g. [8], or strong electron correlations, e.g. 
[9] . Different scenarios refer to different energy scale 
of the charge excitations coupled to BS phonons. 
In most cases its direct experimental measurement 
is not available, because it requires a sufficient mo- 
mentum resolution at relatively low energies. Most 
of the present knowledge comes from optical exper- 
iments, i.e. small momentum q, where the coupling 
to the BS phonon mode is vanishing. In this con- 
text it seems interesting to consider the phonon it- 
self as a probe of the charge excitation. We employ 
a semiphenomenological model where the charge re- 
sponse coupled to the phonon is parametrized in a 
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simple form. To determine all three parameters some 
additional measurements of the BS mode properties 
are required. We propose to consider dependence 
on oxygen isotope substitution of phonon dispersion 
and lincwidth. This would allow not only to deter- 
mine the parameters of the model but also to test 
its validity, since the number of equations is larger 
than required by selfconsistency. 

2. Model and results 

The Dyson equation for the renormalized phonon 
propagator can be written in the form 

luoD- 1 (q, lo)=uj 2 - loI (1 + a sin 2 (q x /2) P (q, u)) , (1) 

where the structure factor corresponds to (100), luo 
is the bare frequency, a is the coupling constant [1] 
and P (q, u>) the charge susceptibility. The latter is 
taken in the Lorentzian form 
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with unknown dispersion Q, qi linewidth T q and os- 
cillator strength rj q . The model (1), (2) defines the 
equation for the phonon frequency 
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and lincwidth 
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where (3 q — ar\ q is the effective coupling constant. 
Analysis of above equations shows that according to 
observed features of the BS softening, dispersion Cl q 
along (100) and (010) should be much smaller than 
in diagonal direction, that agrees with results based 
on t — J model [1] . We then estimate the values of 
T q by substituting known data on u q and 7 9 into (4) 
and changing the values of fl q . After these results 
are confronted with data on optical absorption, it 
follows that the relevant charge excitation should be 
located in the mid infrared region, 0.2 — 0.6 eV. This 
is further corroborated by the analysis of its doping 
dependence in some cuprates [10] as compared to 
softening of the BS phonons. 

Additional information can be extracted by con- 
sidering the effect of the isotope mass shift dm = 
m (O 18 ) - to (O 16 ) on dispersion dui q = uo q (O 18 ) - 



uj q (O 16 ) and linewidth d^ q . In particular, the model 
yields 



dujduo ~ y/l-^sin 2 {q x /2)/^l 
d~f q /duj ~ 2-/ q /ujQ. 
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in the lowest order of the (Vo/£l q expansion. We fur- 
ther define the momentum dependent isotope coef- 
ficient a q = — (mduiq I Lo q dm) that allows to intro- 
duce a new relation to the parameters of the charge 
excitation derived in the same approximation 

a q -a Q = a p q (T 2 q /n 2 q - l) sin 2 (q x /2) uj 2 jn 4 q . (7) 

If is not small, i.e. dynamic effects are impor- 

tant, the full expressions should be used (not given 
here) . 
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Fig. 1. The expected isotope dependence of phonon disper- 
sion along the two directions in the Brillouin zone. 

Eqs. (3-7) contain the relation between experi- 
mentally measurable characteristics of the phonon 
spectrum and the parameters of the model describ- 
ing the charge excitation. On the basis of avail- 
able data the model predicts that the energy of 
charge excitation responsible for the renormaliza- 
tion of the BS phonon mode is located in the mid 
infrared region. Its energy is lowered with hole dop- 
ing and respectively its effect on the BS phonon in- 
creases. To be noted that the sign of the differential 
isotope coefficient (7) depends on the ratio of il q 



2 



to the linewidth. We expect that the fingerprint of 
the broad MIR structure will show up in a positive 
value of this coefficient. Another observation is that 
phonon linewidth could be a more sensitive probe 
of the isotope effect than phonon dispersion as seen 
from (5) and (6). 
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